AD-A265 

IBM 


3lt 


UMENTATIO&  PAGE 


:0  rm  Ippro.e'J 
OMB  Q?Cd  0*8- 


I'*'?  ••-- •*‘J  •  ■■  ir -**  ',*-*-3  4,:  „  f,,  ., 

'"i  3yi*1'  ■.*«*  •’•wQ-.i'-n  i-'.c-i  i'.+  :o#'»r  -r$  **3  a*oc.^\ 

■»  *  :  *3  r  ;c  „-tcc  j 


2  REPORT  OATE 

i  REPORT  type  and  OATES  COVERED  '  “  - 

FINAL/01  SEP  91  TO  31  DEC  92 

4.  TITLE  AWO  SUSTITLf 

"A  SIMPLE  MODEL  OF  MELT  FRACTURE"  (U) 


6  AUTHOR(S) 

Professor  James  Greenberg 


7.  PERFORMING  ORGANIZATION  NAME'S  A  NO  AJOa£5ji£$j 

University  of  Maryland 
Dept  of  Math  &  Statistics 
Baltimore,  MD  21228 


S  FUNDING  NUMBERS 


2 3 04/ A 1 
AF0SR-91-0352 


8  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


ABttR-Tk-  9  3 


9.  SPONSORING /MONITORING  AGENCY  N-V-  «: 
AFOSR/NM 

110  DUNC/aN  AVE,  SUTE  B115 
BOLLING  AFB  DC  20332-0001 


11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

12b.  DISTRIBUTION  CODE 

APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  IS  UNLIMITED 

UL 

13.  ABSTRACT  (Maximum  200  words) 


The  PI  produced  an  excellent  explanation  of  the  unpleasant  shark- skinning  observed 
in  certain  polymer  extrusion  processes.  This  work  has  been  brought  to  the 
attention  of  researchers  at  Corning  and  Hoechst  Celanese  and  Greenberg  and  Demay 
will  work  this  summer  with  members  of  the  Materials  Sciences  Center  at  the  Ecole 
Nationale  Superieure  des  Mines  de  Paris  led  by  J.F.  Agassant.  One  goal  of  this 
work  is  to  see  if  the  same  oscillatory  phenomena  is  present  when  one  replaces  the 
siip  boundary  condition  by  a  no  slip  one  and  looks  instead  at  materials  whose  shear 
stress  -  strain  rate  constitutive  equation  has  a  spinodal  type  nonlinearity.  A 
difficult  question  also  worth  pursuing  is  whether  now  understanding  the  nature  of 
the  flow  instability  -  a  switch  from  a  slip  to  a  no  slip  boundary  condition  at  the 
wall  of  the  capillary  tube  -  if  it  is  possible  to  control  the  inlet  flow  to  the 
capillary  in  the  unstable  regime  in  such  a  way  as  to  reduce  the  oscillations  and 
shark  skinning  of  the  final,  product. 


14.  SUMtcrr  t comic 

IS.  NUMBER  OF  PAGES 

G  !'-- 
O  *-* 

% 

S' 5^ 

>  >  *  u* 

K  t* 

■  *  i  \J 

c 

§ 

* 

IB.  PRICE  CODE 

17.  SECURITY  CLASSIFICATION 

18.  SECURITY  CLASSIFICATION 

19.  SECURITY  CLASSIFICATION 

20.  UftMTATION  OF  ABSTRACT 

OF  REPORT 

OF  THIS  PAGE 

OF  ABSTRACT 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

SAR(SAME  AS  REPORT) 

NSN  7540.01-280-5500  Standard  form  29*  (R#v  2-89) 


»r*w«M*  b*  ANN  tt* 

m-ma  j 


Final  Technical  Report  on  the  Research  based  on  the  proposal  AFOSR- 
91-0352  for  the  period  September  1,  1991  through  December  31,  1992. 

The  original  proposal  was  to  support  two  months  of  Greenberg's 
research  efforts  on  Lattice  Gases  and  Transport  Processes  while  he 
was  on  sabbatical  leave  in  Europe  from  the  period  January  1,  1992 
through  December  31,  1992.  Before  leaving  for  Europe , Greenberg 
contacted  his  program  officer,  Dr.  Arje  Nachman,  and  requested  a 
change  in  the  scope  of  the  project.  Greenberg  proposed  to  work 
instead  on  problems  in  Plastic  Flow  and  Oscillatory  Flows  in  Molten 
Polymers.  Dr.  Nachman  approved  this  request. 

During  the  period  June  1,  1992  through  September  30,  1992 

Greenberg  was  resident  at  the  University  of  Nice  and  worked  with 
Yves  Demay  and  Anne  Nouri  on  the  revised  research  program.  Two 
papers  resulted  from  this  effort: 

(1)  with  Anne  Nouri,  Antiplane  Shearing  Motions 
of  a  Visco-Plastic  Solid;  to  appear  in  SIAM 
J.  Math.  Analysis 

(2)  with  Y.  Demay,  A  Simple  Model  of  Melt  Fracture,  to 
appear  in  European  Journal  of  Applied  Mathematics. 

The  latter  paper  nicely  compliments  work  of  Greg  Forest,  also 
supported  in  Dr.  Nachman's  program,  and  gives  an  excellent 
explanation  of  the  unpleasant  shark-skinning  observed  in  certain 
polymer  extrusion  processes.  This  work  has  been  brought  to  the 
attention  of  researchers  at  Corning  and  Hoechst  Celanese  and 
Greenberg  and  Demay  will  work  this  summer  with  members  of  the 
Materials  Sciences  Center  at  the  Ecole  Nationale  Superieure  des 
Mines  de  Paris  led  by  J.  F.  Agassant.  One  goal  of  this  work  is  to 
see  if  the  same  oscillatory  phenomena  is  present  when  one  replaces 
the  slip  boundary  condition  by  a  no  slip  one  and  looks  instead  at 
materials  whose  shear  stress  -  strain  rate  constitutive  equation 
has  a  spinodal  type  nonlinearity.  A  difficult  question  also  worth 
pursuing  is  whether  now  understanding  the  nature  of  the  flow 
instability  -  a  switch  from  a  slip  to  a  no  slip  boundary  condition 
at  the  wall  of  the  capillary  tube  -  if  it  is  possible  to  control 
the  inlet  flow  to  the  capillary  in  the  unstable  regime  in  such  a 
way  as  to  reduce  the  oscillations  and  shark  skinning  of  the  final 
product . 


c 


DTIC  QTJALJTy  INSPECTED  Si 


A  Simple  Model  of  Melt  Fracture 


J.  M.  Greenberg* 

Department  of  Mathematics  and  Statistics 
University  of  Maryland,  Baltimore  County 
Baltimore,  MD  21228 

Y.  Demay 

Institut  Non-Lineaire  de  Nice 
U.M.R.-C.N.R.S.  nol9 
Universite  de  Nice 
Parc  Valrose,  06034  Nice  (France) 


1.  Model  Development 


Accesion  For 

NTIS  CRA&l 
DTIC  TAB 
Unannounced  □ 
Justification 


By . . 

Distribution/ 


Availability  Codes 


Dist 

Avail  and/or 

Special 

M 

A  classical  measure  of  viscosity  of  a  molten  polymer  is  obtained  by  studying  the  re¬ 
lationship  between  the  flow  rate  and  pressure  drop  when  the  polymer  is  flowing  in  a 
capillary  tube  of  radius  R  and  length  L.  For  a  linearly  viscous  polymer  with  viscosity 
Tjj  this  relationship  takes  the  form: 


(1.1) 


R2  A  P 
Srji  L 


where  A P  is  the  pressure  drop  across  the  capillary  from  z  =  0  to  z  —  L  and  w  is  the 
mean  velocity  taken  over  a  cross-section  perpendicular  to  the  flow  direction.  For  such 
a  flow,  the  wall  shear  is  given  by 


(1.2) 


RAP 

2L 
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Figure  1. 


and  thus  (1.1)  is  equivalent  to  the  following  relationship  between  w  and  tw: 

(1.3)  w  =  — — tw. 

4??i 

These  relationships  follow  from  the  local  equilibrium  equations  for  a  linearly  viscous 
fluid  subjected  to  a  constant  pressure  gradient  and  a  no  slip  boundary  condition  at  the 
capillary  wall.  The  underlying  flow  in  this  case  is  Poiseuille.  If  the  no  slip  boundary 
condition  is  replaced  by  the  slip  boundary  condition 


(1.4) 


w(R)  = 


R(4F  —  1) 

4?7i 


FSi' 


then  the  relation  (1.2)  still  obtains  but  (1.1)  is  replaced  by 


(1.5) 

or  equivalently 

(1.6) 


w 


w 


R2FAP 

2r/il 

RF  , 


7i 


^hese  relations  are  derived  later  in  this  section;  for  details  see  ( 1 .27)— ( 1 .32). 
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Pressure  Drop 
D* 


Mean  Flow 
Figure  2. 


For  molten  polymers  one  typically  observes  two  distinctly  different  stable  steady 
flow  regimes  (see  Figure  2).  For  mean  velocities  w  <  Wi  one  operates  on  the  curve 
5i  whereas  for  mean  velocities  w  >  ui2  one  operates  on  the  curve  £2.  Experimental 
evidence  indicates  that  the  pressure  drop  D\  associated  with  the  velocity  wx  is  larger 
than  the  drop  D2  associated  with  w2.  In  general  one  finds  no  stable  steady  flows  with 
velocities  between  Wi  and  w2;  rather  one  observes  oscillatory  flows  with  mean  velocities 
in  the  range  (u>i,  w2)  and  pressure  drops  in  the  range  (D2,  D\ ),  One  normally  interprets 
such  data  as  follows:  on  the  curve  S\  the  polymer  satisfies  a  no  slip  boundary  condition 
at  the  capillary  wall  whereas  on  the  curve  S2  slipping  at  the  wall  is  taking  place.  In 
general,  the  curves  S\  and  S2  are  nonlinear  which  indicates  that  the  polymer  is  not  a 
simple  linearly  viscous  fluid  and  that  the  slip  boundary  condition  is  really  modeled  by 
a  nonlinear  relationship  between  w(R )  and  tw.  These  details  will  not  concern  us  here. 
Our  goal  is  a  simple  model  which  explains  the  gross  features  of  the  transition  from 
steady  to  oscillatory  flows. 

Motivated  by  (1.3)  and  (1.6)  we  assume  that  w  and  rw  are  related  by 


(1.7) 


w  = 


RF 


where  rji  is  the  viscosity  of  the  polymer,  R  is  the  radius  of  the  capillary,  and  F  is  a 
dimensionless  parameter  which  ranges  between  1/4  and  F2  >  1/4.  Bearing  in  mind  the 
relationship  (1.2),  we  assume  the  existence  of  a  switch  curve  (scaled  as  shown) 


D\  R 
2  L  ’ 


IV  < 


(1.8)  S(w)  = 


£i*  +  &( . , ..Pi- Pi, . ^  (,7,  _  #DA 

2 L  ~  2L  [  R^FtDi  rTdT  I  \w  Svi  L  )  ' 

\  L  L  / 


R2Di 

8r?|6 


DiR 
2  L  1 


r 
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which  determines  the  stability  of  the  curves  S\  and  S2\  specifically  we  assume  that  F 
switches  from  1/4  to  F2  via  the  transition  rule 

nox  dF(F2-F;  |ru|>%)  =  5(FM 

1  Jo  dt  ~  \  }  -  K  W  <  S(w)  =  S(F\tw\) 

where  To  is  a  fixed  relaxation  time  (see  Figure  3). 


Equation  (1.7)-(1.9)  are  constitutive  equations  relating  u>,  F,  and  and  these  hold 
at  each  point  z  along  the  capillary.  To  these  we  must  adjoin  variants  of  the  continuity 
and  balance  of  momentum  equations  together  with  an  equation  of  state  relating  the 
density  to  the  pressure  in  the  polymer.  Again,  the  geometry  of  our  system  is  shown 
in  Figure  1  and  our  model  accounts  for  what  happens  in  the  capillary  tube,  not  the 
reservoir. 

We  assume  the  flow  in  the  capillary  is  axisyrnmetric  and  that  all  quantities  depend 
only  on  r,  z,  and  t.  We  let  p  denote  the  polymer  density,  p  denote  the  pressure,  and 
assume  that  the  velocity  field  u  and  Cauchy  stress  tensor  a  are  of  the  form: 

(1.10)  u  =  uer  +  we^ 

and 

(1.11)  a  -  — p(er  ©  er  4-  e8  ©  e9  +  e#  ©  £3)  +  avTrtT  ©  er  +  cr99e9  ©  e9 

+  <?33&J  ©  §3  +  ^(£3  ©  er  +  er  ©  63). 2 

2er  =  (cos0,sin0,O),  S.g  —  (—  sin  cos  0,0),  and  for  any  vectors  a  =  (01,02,  03)  and  b  =  (61,621^3) 

(  \ 

we  let  a  ©  b  —  aTb  where  aT  —  I  02  1 . 

\  a3  / 


f 
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(a?r,  (Tg3,  cr$. j,  r)  are  the  nonzero  components  of  the  viscous  stress  tensor  and  are  assumed 
related  to  u  and  tv  by 


(1.12) 


rv  — 


A,(if(ru)  +  |S)  +  2„& 

A,(l|M  +  |)  +  2„  = 

^  =  A,(i|(rt)  +  te)  +  2„t 
^  +  £)• 


u 


We  do  not  make  the  Stokean  hypothesis  3Aj  +  2771  =  0  but  we  do  assume  that  Aj  = 
<9(t7i).  The  governing  equations  are  the  continuity  equation  and  balance  of  momentum 
in  the  directions  er,  e^,  and  e 3.  These  take  the  form 


(C) 


dp  1  d 
dt  ^  r  dr 


(rpu)  +  pw ^ 


=  0, 


(M,) 

and 

(M3) 


1  d 


(Mr)  -Qt(pu)  +  zi*z(rPu  )  +  itzipu™)  ~  ~zz  +  + 


dp  ,  d<Tvrt  «  -  alg)  dr 


dr 


dz 


dr  dr 

_La?  =  0 

r  d9  ’ 


+  Tz~ °’ 


1  d 


d 


dp  ,  dav32  1  d 


At  the  boundary,  r  —  R,  we  assume  that 


(BC) 


u(R,z,t)  =  0  and  z,£) 


4t/i 


(4F-  l)rw(z,f), 


where  F  >  1/4  is  a  dimensionless  quantity  depending  on  z  and  t  and  rw(z,  t )  =  r(R,  z,  t) 
is  the  wall  shear.  In  order  that  the  stresses  <7/,.,  crgB,  and  3  be  finite  at  r  =  0,  we  also 
require  that  the  radial  component  of  the  flow,  u,  vanishes  faster  than  r  as  r  tends  to 
zero.  To  close  the  system  we  assume  that  the  polymer  is  slightly  compressible  and  that 
the  following  equation  of  state  holds  between  the  density  and  pressure: 


(EOS) 


P 


Po 


Here,  pQ  and  p0  are  reference  values  of  the  density  and  pressure  and  0  <  t\  <C  1  is  a 
dimensionless  small  parameter. 

To  assess  which  terms  in  this  system  are  important  and  which  may  be  neglected  we 
cast  the  system  in  dimensionless  form.  We  let 


r  =  Rr  1,  z  =  Lz\, 


and 


_  2et7i£2 

poR2 


(1.13) 


f 


auu 

P  =  PoPi ,  W  =  u  =  jjjfjiii,  p  =  2popi, 

(1-14)  -  r  =  2af ->  <r  =Po(f)2<r,>  ^=Po(f)2^fl, 

.  ff33  =  i*  (l)2  *33»  and  T«>  =  EiTrJ,- 

Then  it  is  a  relatively  easy  calculation  to  show  that  (1.12)  transforms  to 

"Ii.i  =  /<  (nsr(r‘»i) + 1?) +  2I? 

,11M  +  fc)+2?f 

.  r1  =  fe  +  <*) 


(1.15) 


where 

(1.16) 


f,  =  —  =  0(1)  and  e2  =  t  <  1 

Pi  L 


and  the  equation  of  state  becomes 

(1.17)  pi  =  (1  +  2eipi). 

The  continuity  equation  transforms  to 


U-18) 


1  9pi  1  9  ,  ,  d  ,  .  n 

2irar+W,'”u,)+^('’,“")=0 


and  (1.18),  when  combined  with  (1.17),  yields 


|,+1a(p)  +  M+2£i 

9<i  rx  or  i  ozi 


(  1  ^(riPiMi)  d(piwi) 

V !  9rq  <9?! 


Finally  the  three  balance  laws  (Mr),  (Ms),  and  (M3)  take  the  form 

m  om  4*3 d(pim)  ,  4t 3  5  , _ 2\  ,  .4.  d(pl«l»l) 

(1-20)  2^-dir  +  —d^iiPiU')+V3~d^— 


(1.21) 


«r,  -  <4?)  ,  dor' 


I^i  =  o 

r,  90  ’ 


,  ,2^1!  _  A 

+  2  dzi  2  dr  1 ' 


(1.22) 


f|f3  djpiWi) 

2et  dti 

o  9^  o. 


2  /  1  9(r,piU!iui)  9(piWj)\ 

+  e’f3  In'  Sr,  +  -  &T j 

o.  _  1  1  d{riTl) 


c4(l  -F  2eipi)  =  2- 


7 


where  ex  is  the  small  parameter  in  (EOS),  e2  =  jr  <.  1  and 

j  /  /)a  TVi  /?2 

(1.23) 


def  PoPO#  ,  POST 

e3  =  - , —  and  e4  = 


•7? 


Po 


2 

Our  principal  simplification  comes  from  assuming  that  and  £4  are  small  and  from 
neglecting  all  terms  with  c’s  in  ( 1 . 19)— ( 1 .22).  This  leads  us  to  the  following  reduced 
system 

.  .  dp  1  1  d(rxux)  dw\ 

(1.24  ~  + - U-1Z  +  — J=0, 

dtx  rj  drx  dz  1 


(1.25) 

and 

(1.26) 


dpi  _  dpi 

drx  dO 


=  0, 


1  ^(rir1)  _  ^dpx 
rx  drx  dzx 


Equation  (1.25)  implies  that  px  =  pi(zx,tx)  and  thus  (1.26)  yields 

(1.27) 

From  the  last  relation  we  see  that  the  scaled  wall  shear,  —  Ti(l,Zi,fi),  is  given  by 

(1.28) 


1  dPi 

r  =  rj- — . 
OZ\ 


dpi 


dzx 

The  constitutive  equation  for  r1  when  t2  =  0+: 

dv)X 
dr. 


(1.29) 

together  with  the  scaled  form  of  (BC)2: 

(1.30) 


r1  = 


„  ,  (4F-1)  , 

u>i(l,zi,ti)  = - - - rw 


implies  that 

(1.31) 

and  that 

(1.32) 


wX{ri,zi,tx)  =  Th  0  <  n  <  1, 


wi{zi,ti)  =  2  f  wx{rx,zx,tx)rxdrx  =  -Fr*. 
Jo 
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If  we  multiply  (1.24)  by  iq,  integrate  the  resulting  equation  from  r,  =  0  to  r,  =  1, 
and  exploit  the  fact  that  rtui  vanishes  at  rt  =  0  and  rj  =  1  we  obtain  the  following 
equation  for  p\{z,tx)  and  tD(ij,<i): 


(1.33) 


dpi  dwx 

dtx  dzx 


From  this  Ictst  identity  and  (1.24)  we  also  obtain 


(1.34) 


1JL 

ri  drx 


(n«i) 


d 

dzx 


(wj  -  Wi) 


where  wx  and  wx  are  given  by  (1.31)  and  (1.32).  Equation  (1.34)  may  be  integrated  by 
quadrature  yielding 


(1.35) 


d  /I  fTl~ 

Ul(ri,2i,ti)  =  —  (  —  /  (wX(zX,tX)~  WX(s,Z  1,<J 
UZ\  \  V |  Jo 


The  above  equation  guarantees  that  as  rq  — ►  0+,  — ' tl- converges  to  5577(^1(^1,  <i)~ 
u?(0,  zj,<)),  a  finite  limit.  The  existence  of  such  a  limit  was  required  in  order  that  the 
stresses  c-”r,  av&6,  and  were  finite  at  the  center  line  of  the  capillary.  We  also  note  that 
fj(w(zx,tx)  —  wx(s,  zi,  tx))s  ds  —  0  and  thus,  as  defined,  does  meet  the  boundary 
condition  linv,.,!- ux(rx,  zx,  tx)  =  0. 

To  summarize,  our  remaining  equations  for  px  and  wx  are  (1.28),  (1-32),  and  (1.33). 

To  these,  we  adjoin  the  scaled  versions  of  (1.8)  and  (1.9)  for  the  evolution  of  F.  We 
write  the  constants  Dx  and  D2  of  (1.8)  as 


(1.36) 


Dx  =  Pqcx  and  D2  =  pQc2 


and  observe  that  (1.8)  and  (1.9)  transform  to 


(1.37) 


xdF_tF2-F;  |ri|  >$(*,) 
dtx  |t^,|  <  S(w\) 


where 

(1.38) 


S{wx)  =  <; 


Cl, 

Cl  + 
C2, 


(ca  — ci ) 
(fscj-ci/4) 


(WX  -  Ci/4), 


wx  <  Ci/4 

^  <wx  <  F2c2 
F2c2  <  wx 


and 

(1.39) 


T0R2po 

expxL2 


and  again  T0  is  the  fixed  relaxation  time  appearing  in 


(1.9).  These  equations  may  be 
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combined  to  yield  the  following  system  for  px  and  F 


(D) 

fa  _  A  fr^l)  -  n 
dt,  dz,  (  dz,) 

and 

dF  {F-F  -fe>5(-Ffe) 
ai‘  H-F  -^<S{-F^) 

(F) 

where  again  S  is  given  by  (1.38).  This  latter  system  is  solved  in  0  <  zx  <  1,  the 
normalized  region  of  the  center  line  of  the  capillary,  subject  to  the  following  boundary 
conditions: 

(BC)  -  F(0,  =  q  >  0  and  pi(Ui)  =  1- 

Oz\ 


The  condition  at  zx  =  0  corresponds  to  the  fact  that  the  input  of  material  from  the 
reservoir  to  the  pipe  is  done  at  a  constant  rate  and  the  latter  condition  expresses  the 
fact  that  at  the  pipe  exit  the  material  is  at  atmospheric  pressure. 


2.  Analysis  of  (D),  (F),  and  (BC) 


Our  interest  here  is  in  the  analysis  of  the  system 


(D) 


dp 

dt 


0  <  z  <  1 


(F) 


FW.  -£>F(-F§?) 
i-F,  -ff  <  S  (-F|f) 
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where  again 

(S) 

S(  w)  =  < 

f  Cl, 

1  Cl  4.  .  a) 

I  1  ^  (F2c2-c,/4)'.U; 

l  c2, 

w  <  ci/4 
ci/4  <  w  <  F2c2 
F2c2  <  u 

These  equations  are  solved  subject  to  the  boundary  and  initial  conditions 
(BC)  -  F(0,t)|?(0,i)  =  q  >  0  and  p(l,<)  =  1 

0Z 

and 

(IC)  p(z,0+)  =  Po(z)  >  0  and  F(z,0+)  =  F0(z)  €  (1/4,  F2),  0  <  z  <  1. 

In  this  section  we  drop  the  subscript  1  for  normalized  dimensionless  quantities. 

It  is  not  particularly  difficult  to  demonstrate,  either  analytically  or  computationally, 
that  if  the  parameter  q  in  (3C)  satisfies  q  <  c{/4,  then  solution  of  (D),  (F),  (BC),  and 
(IC)  converges  to  the  stable  equilibrium 

(2.1)  Poo(z)  =  1  +  49(1  -  z)  and  F^(z)  =  1/4,  0  <  z  <  1, 

whereas  if  q  satisfies  q  >  F2c2 ,  then  p  converges  to  1  +  ?(1  —  z)/F2  and  F  to  F2.  Thus 
we  shall  confine  our  attention  to  the  situation  where  q  €  (ci/4,F2c2)  for  it  is  in  this 
range  that  we  see  the  oscillatory  flows. 

We  use  an  operator  splitting  technique  to  integrate  the  system  (D),  (F),  (BC),  and 
(IC).  If  6  is  our  time  step,  we  assume  we  have  an  approximate  solution  (p",Fn)(z) 
defined  at  tn  =  nS,  n  =  0,1,....  During  the  first  half  step  we  advance  p  keeping  F 
fixed;  that  is,  we  solve 


MW)-- 


subject  to  the  boundary  condition  (BC)  and  the  initial  condition 


(2.3)  (p,F)(z,0)  =  (p",Fn)(z),  0  <  z  <  1. 

We  denote  the  solution  to  (2.2),  (BC),  and  (2.3)  at  t  =  8  by  (p1^2,  Fl^2)(z)  and  of 
course  F’^2(z)  =  Fn(z).  During  the  second  half  step  we  solve 

al_n  ,  .acfn-F,  -|f  >s(-f§‘) 

st  »  H-f.  -sf<5(-e|f) 

subject  to  the  initial  condition 

(p,F)(z,0)  =  (pl'!,F")(z),  0  <  z  <  1. 


(2.5) 
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The  solution  to  this  latter  problem  at  t  =  6  is  the  new  approximate  solution  at  time 
(n  +  1)5. 

In  solving  (2.2)  and  (2.4)  we  let  n  be  an  integer,  h  —  5^7,  and  evaluate  p  and  F 
at  the  respective  grid  points  z\  —  ^~h  and  z\  =  (k  —  1  )h\  1  <  k  <  (n  +  1).  We  use 
a  fully  implicit  first  order  method  on  both  (2.2)  and  (2.4).  This  allows  us  to  choose  6 
and  k  independently  and  have  a  stable  integration  scheme. 

We  have  run  a  number  of  simulations  on  this  system  with  a  variety  of  parameter 
values  (Ft,  <70,  C2,  Cj,  A)  and  shall  report  on  a  few  of  these  at  the  end  of  this  section. 
Before  doing  this  though,  we  report  on  one  finding  which  surprised  us.  We  observed 
that  once  regular  oscillations  had  been  established  (and  transients  had  died  out)  that 
the  pressure  field  was  approximately  linear  at  each  instant  of  time;  that  is  the  p  profile 
was  approximately  given  by 

(2.6)  p  ~  pun  po(<)  +  (1  -  po(0)«- 


This  observation  suggested  that  a  simplified  model  of  the  oscillatory  flows  should  be 
possible.  Below  we  present  such  a  model  which  exploits  (2.6). 

If  we  integrate  (D)  from  2  =  0  to  2  =  1  and  exploit  the  fact  that  the  solution  to 
(D)  is  approximately  given  by  (2.6),  we  find  that  po  must  satisfy 


(2.7) 


1  dp0 

2  dt 


—  (F(l,t)  —  F(0,t))(\  —  po(t)). 


dp  dpn 

Here,  we  are  exploiting  —  ~  "  =  (1  —  po).  A  consistent  interpretation  of  the 

,  ,  .  az  dz 

boundary  condition 


(BC) 

is  that 

(2.8) 

and  thus  (2.7)  reduces  to 
(2.9) 


F(0,f)g(0,<)  =  ? 


-  F(0,t)(l  -  po(t))  =  q 

1  dpo 

2  dt 


=  q~  F(\1t)(po{t)  -  1). 


dP, 


The  equation  for  F(l,t)  is  obtained  from  equation  (F);  we  merely  replace  —  ^(1,<)  by 

oz 


^Plin 

dz 

(2.10) 


(1,<)  =  (po(0  —  1)  and  obtain 


F2  -  F(l,t),  (po(t)  -  1)  >  5(F(1,0(po(0  -  1)) 

4  ~  F(l,<),  (po(t)-l)<S(F(l,t)(po(t)~l)). 
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The  system  (2.9)  and  (2.10)  is  the  simplified  model  used  to  describe  the  oscilla¬ 
tory  flows.  For  each  q  £  (^-,c2F2),  the  solutions  of  (2.9)  and  (2.10)  converge  to  a 
unique  orbitally  stable  limit  cycle  whose  particular  character  depends  on  the  choice  of 
parameters  (\,q,  F2,ci,c2).  Below  we  show  two  figures.  Both  were  run  with 

(2.11)  (A1,F2,c1,c2)  =  (.5,.3,L95,1.85) 

In  Figure  5 

(2.12)  q  =.7 
and  in  Figure  6 

(2.13)  q  =  .6. 

In  the  top  portrait  of  each  Figure  the  horizontal  axis  represents  the  exit  flow  velocity, 
Wexit{t),  which  is  computed  by 

(2.14)  W'xu(t)  =  ~  1) 

and  the  vertical  axis  represents  the  pressure  drop  across  the  tube,  po{t)  —  1.  The  closed 
curve  in  each  figure  is  the  stable  limit  cycle  for  the  indicated  parameter  values  and 
these  demonstrate  the  sensitivity  of  the  oscillations  to  the  incoming  value  of  the  flow, 
q.  The  bottom  graphs  of  each  Figure  show  the  time  histories  of  the  pressure  drop  and 
exit  velocity. 

Figure  7  represents  a  simulation  of  the  full  system  (D),  (F),  (BC),  and  (IC)  run  with 
the  same  parameter  values  used  to  generate  Figure  6.  In  these  simulations  n  =  100 
(and  thus  dx  ~  .01)  and  dt  =  .005.  In  the  top  portrait  of  Figure  7  we  have  shown  (1) 
the  limit  cycle  for  the  reduced  system  with  the  parameter  values  used  for  Figure  6,  (2)  a 
graph  of  the  pressure  drop  (p(0,  t)~  1)  versus  the  exit  velocity  u?(l ,  <)  =  —  F(l,  i)|£(l,  t), 
and  (3)  a  graph  of  minus  the  exit  pressure  gradient,  — ||(l,t),  versus  the  exit  velocity 
=  —  F(l, f)§|(l,  t).  The  pressure  drop  curve  is  the  narrower  of  these  two.  The 
second  set  of  pictures  in  Figure  7  shows  time  histories  of  these  quantities.  In  the  graph 
labeled  “pressure  drop  vs  time”  the  time  history  of  —  ||(l,t)  is  the  curve  with  the 
largest  oscillations. 

Figures  8-39  show  snapshots  of  relevant  flow  variables  at  the  time  indicated  in  each 
figure.  The  data  shown  in  Figure  8  was  also  used  to  generate  Figure  7.  The  only  pictures 
requiring  some  explanation  are  the  ones  labeled  “shear  vs  velocity”.  Recalling  that  the 
normalized  wall  shear,  tw  <  0,  satisfies  =  — | \{z,t)  and  that  the  normalized  mean 
velocity,  u?(z,f),  is  given  by  u>(z,t )  =  —F(z,t)&(z,t),  we  obtain  the  graphs  under 
discussion  by  plotting  the  curves  2  — ►  (  — F(z,  Of^z,  <),  — §^(z,  f)),  0  <  z  <  1,  at  the 
times  indicated.  The  points  labeled  x  are  the  images  of  z  =  1  and  those  labeled  o  are 
the  images  of  z  =  0  at  the  indicated  times. 
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3.  Conclusions 

We  have  produced  a  robust  model  which  is  capable  of  producing  relaxation  oscillations 
which  are  qualitatively  similar  to  those  observed  in  polymer  extrusion  experiments. 
The  heart  of  the  model  is  the  switch  rule  (1.8)  and  (1.9)  which  allows  us  to  go  from  the 
slip  to  the  no  slip  boundary  condition  at  the  capillary  wall.  Though  we  have  not  done 
it  here,  we  believe  it  is  possible  to  go  from  the  linear  constitutive  equations  employed 
here  to  more  realistic  ones  for  a  molten  polymer  and  to  modify  the  linear  slip  boundary 
condition  to  one  employing  a  nonlinear  relation  between  the  velocity  at  the  wall  and 
the  the  wall  shear  which  is  scaled  by  a  dimensionless  factor  F  which  obeys  a  transition 
mechanism  similar  to  (1.8)  and  (1.9).  Our  colleagues  at  the  Ecole  Nationale  Superieure 
des  Mines  de  Paris  are  currently  pursuing  this  work. 
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ANTIPLANE  SHEARING  MOTIONS 
OF  A  VISCO-PLASTIC  SOLID* 

J.  M.  GREENBERG'  AND  ANNE  NOL'RI' 

Abstract.  The  authors  consider  antiplane  shearing  motions  of  an  incompressible  isotropic 
visco-plastic  solid.  The  Sow  rule  employed  is  a  properly  invariant  generalization  of  Coulomb  sliding 
friction  and  assumes  a  constant  yield  stress  or  threshold  above  which  plastic  flow  occurs,  hi  tins 
model  stresses  above  yield  are  possible;  but  when  this  condition  obtains,  the  plastic  flow  rule  forces 
the  plastic  strain  to  change  so  as  to  tower  the  stress  levels  in  the  material  and  dissipate  energy.  On 
the  yield  surface,  the  flow  rule  looks  like  the  classical  one  for  a  rate  independent  elastic-perfectiy 
plastic  material  when  the  velocity  gradients  are  small  enough  but  differs  from  the  classical  model  for 
large  gradients. 

Key  words,  plastic  waves,  visco-plasticity,  time-dependent  problems 

AMS  subject  classifications.  73E70.  73E6Q.  73E50 


1.  Introduction.  In  this  note  we  consider  antiplane  shearing  motions  of  an  in¬ 
compressible  isotropic  visco-plastic  solid.  This  work  generalizes  and  compliments 
earlier  work  of  Greenberg  [I],  [2],  where  he  considered  simple  shearing  flows  for  such 
materials.  The  flow  rule  we  employ  is  a  properly  invariant  generalization  of  Coulomb 
sliding  friction  and  assumes  a  constant  yield  stress  or  threshold  above  which  plastic 
flow  occurs.  As  with  most  such  theories,  we  assume  a  multiplicative  decomposition  of 
the  deformation  gradient  into  an  elastic  and  plastic  part,  and  we  assume  further  that 
the  deviatoric  part  of  the  Cauchy  Stress  tensor  depends  only  on  the  elastic  portion  of 
the  deformation  gradient.  For  antiplane  shearing  motions  this  decomposition  presents 
no  precedence  problems;  i.e.,  does  the  elastic  deformation  precede  the  plastic  or  vice 
versa?  One  key  feature  of  this  model  is  that  stresses  above  yield  are  possible.  When 
this  condition  obtains,  the  plastic  flow  rule  forces  the  plastic  strain  to  change  so  as 
to  lower  the  stress  levels  in  the  material  and  dissipate  energy.  The  principal  difficulty 
in  formulating  this  model  occurs  when  the  stress  is  at  yield.  Motivated  by  results  of 
Seidmaa  [3],  Utkin  [4],  and  Filippov  [5]  on  sliding  modes  induced  by  discontinuous 
vector  fields,  we  are  led  to  the  flow  rule  advanced  in  (2.38).  On  the  yield  surface,  this 
flow  rule  looks  like  the  classical  one  for  a  rate  independent  elastic-perfectiy  plastic 
material  when  the  velocity  gradients  are  small  enough  but  differs  from  the  classical 
model  for  large  gradients.  This  rule  differentiates  between  loading  and  unloading  and 
generates  an  energy  identity  which  guarantees  that  uniqueness  obtains  for  initial  and 
initial-boundary  value  problems. 

The  organization  of  this  paper  is  as  follows.  In  §2  we  develop  the  appropriate 
equations  describing  antiplane  shearing  flows  in  visco-plastic  solids.  Section  3  focuses 
on  the  uniqueness  issue.  Our  basic  estimate  is  that  the  energy  associated  with  the 
difference  between  two  solutions  generated  by  the  same  data  is  nonincr^asing.  This 
estimate  relies  in  an  essential  way  on  the  definition  of  the  plastic  flow  rule.  In  §4  we 
examine  a  one-dimensional  signalling  problem  and  discuss  (1)  the  structure  of  this 
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solution,  and  (2)  a  procedure  to  analytically  obtain  an  approximate  solution.  We 
also  compare  this  solution  with  what  obtains  for  the  more  studied  model  of  a  rate 
independent  elastic-perfectlv  plastic  material  where  uniqueness  fails.  Section  5  deals 
with  a  numerical  experiment  for  a  two-dimensional  signalling  problem  in  the  corner 
domain  r  >  0  and  tr/2  <  9  <  2sr.  Here  the  stresses  are  singular  as  one  approaches  the 
corner  and  care  must  be  taken  in  the  implementation  of  the  boundary  conditions. 

We  note  that  in  the  last  several  years  there  have  been  a  number  of  other  efforts 
aimed  at  capturing  the  essence  of  plastic  flows.  Antman  and  Szymczak  [6],  [7]  have 
advanced  a  finite  deformation  theory  of  such  materials  which  is  similar  in  spirit  to  ours 
but  differs  in  a  number  of  essential  ways.  Their  model  is  formally  rate  independent 
where  ours  is  not  but  their  model  also  requires  a  history  dependent  strain  hardening 
mechanism.  The  predictions  of  the  two  theories  are  often  qualitatively  different;  these 
differences  arise  since  in  their  model  the  imposition  of  large  loads  tends  to  elevate  the 
yield  stress  and  create  a  temporally  constant  permanent  plastic  deformation,  whereas 
in  our  model  such  loading  would  generate  a  constant  plastic  deformation  rate  and  thus 
a  plastic  deformation  which  varies  linearly  in  time.  This  may  be  seen  by  examining 
the  solution  constructed  in  §4.  Other  efforts  on  elasto-plastic  modelling  may  be  found 
in  Coleman  and  Owen  [8],  Buhite  and  Owen  [9],  Coleman  and  Hodgdon  [10],  and 
Owen  [11]. 

2.  Model  development.  We  say  that  a  body  is  undergoing  antiplane  shear  if 
material  points  £  =  £iej  +  Ge2  +  foes  move  to  x  =  xiei  +  12^2  +  £3^3  with 

(2-1)  *t=£i,  *2=6.  and  x3  =  6  +  ,&>,  0 

under  the  action  of  a  Cauchy  stress  tensor  of  the  form 

T  =  — 3-(ei  ®  ei  +  ez  ®  e2  +  e3  <$>  e3) 

(2.2)1  +(S'uei  <85  ej  -f-  Sojej  ®  ej  +  63363  ®  e3) 

+6'3i(e3  ®  e3  +  e3  eg  ei)  +  632(62  ®  e3  +  e3  c*>  e2). 

Here,  x  is  the  hydrostatic  pressure  and  S  is  the  deviatoric  stress  tensor  and  satisfies 


(2-3) 


trace  (6')  =  6'n  +  S22  +  633  =  0. 


Relative  to  the  above  basis,  the  matrix  representation  of  the  Cauchy  stress  is  given 
by 


(2.4) 


•6’n  0  631 

0  622  632 

631  632  633 


and  relative  to  the  same  basis  the  deformation  gradient  is  given  by 
(2.5)  T  = 


1 


1  e,  =  0 


e2 


°  \  0  \ 

1  I  ,  and  e3  =  I  0  1  are  the  standard  basis  elements  for  R 3  and 

\°/  W  \  ■  ) 

ej  Q  ej  as  ejej  T  are  the  standard  basis  elements  for  linear  operators  from  R3  to  R3 . 
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where 

do  do 

(2.6)  FM  =  £  and 

Noting  that  matrices 

/  l  0  0  \ 

d=  (  o  i  o  J 

satisfy  the  commutation  relation 

(2*  0  ,bi)  -^(ai ,6i +67 ) i 

we  feei  justified  in  decomposing  the  deformation  gradient  T  into  its  elastic  and  plastic 
parts  £  and  V  by 


£  :d=f  I 

f  1 

0 

0  \ 

.  .  /  1 

0 

0  \ 

(2.8) 

0 

1 

0 

and  V  :  =  0 

1 

0 

\ 

^  e3\ 

«32 

1  ) 

\  P31 

P32 

1  ) 

where 

(l  0  0\ 

(2.9)  T  =  £V  =  V£  =  0  1  0  J  . 

\  «31  +  P31  632  +  P32  1  / 

For  such  antiplane  shear  flows  one  need  not  make  any  assumption  about  the  prece¬ 
dence  of  the  elastic  and  plastic  parts  of  the  flow. 

Our  basic  constitutive  assumption  is  that  under  a  change  of  reference  frame  E 
transforms  in  the  same  way  as  F  and  that  the  deviatoric  stress  S  is  an  isotropic,  frame 
indifferent,  trace  free  function  of  the  elastic  deformation  gradient  E. 2  The  constraint 
that  5  is  an  isotropic,  frame  indifferent  function  of  E  implies  that  5  must  have  the 
functional  form 

(2.10)  5  =  «/  +  0££T  +  y£~r£~ 1 

or 

(2. 1 1) 3 

/  1  0  0  V  /  1  0  e3l  \  /  I-f«3!  C31«32  ~ «31  \ 

S  =  a|  0  1  0  J  +  3  jo  l  e32  I  +7  I  e3l^32  H-C§2  —  e32  1 

\  0  0  1  )  \  e31  «32  l+e31+e32  /  \  — e31  ~e32  1  J 

where  «.  3,  and  7  are  functions  of  the  invariants  of  ££T ,  in  this  case  the  scalar 
e3i  +  ?§,.  Equation  (2.4)  implies  that  S?  1  =  S12  =  0  and  this,  in  turn,  implies  that 
7  =  0  while  the  condition  that  trace,?  =  0  implies  that  o  =  — ,d(l  +  ((e^,  -F  eln}/.1))). 
Combining  these  identities  with  (2.11)  yields 

/  3(^31  +  e32 )  0 

(2.12)  5  =  /?  0  -^5i+4>) 

\  «31  «32 


C31  \ 

«32  I 

§(e31  +  e32)  / 


3  E  is  the  tensor  whose  matrix  representation  relative  to  the  basis  elements  ej  0  e:  is  given  by 
(2.3),. 

3  For  details  see  Gurtin  [12]. 
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In  the  sequei  we  shall  assume  that  0  is  a  positive  constant.  Equation  (2.12)  implies 
that  we  may  regard  the  elements  S31  and  S32  as  basic  descriptors  of  our  system.  In 
terms  of  these  5  and  £  take  the  form 


(2.13)  5  = 


-^(sh  +  sl,) 


-MSI  +  S& 


(2.14) 


1  0  0 

0  1  0 

S3,  $52  . 

0  0 


We  now  turn  to  the  equations  of  motion.  Equation  (2.1)  implies  that  the  Eulerian 
velocity  field  u  is  of  the  form 


(2.15) 
where 

(2.16) 


u  =  u(zl,r2,<i)e3, 


u(xl,x2,ti)  =  ^-(zi,z2,<i); 


and  (2.16),  when  combined  with  (2.6),  implies  that 


(2.17) 

df3 1 

dtx 

1! 

0 

and 

(2.18) 

dFz2 

dt\ 

0 

II 

<§11 

1 

Additionally,  (2.9)  and  (2.14)  imply  that 

(2.19)  F3l  =  ^-  +  P3, 
and 

(2.20)  F32  =  +  P32- 

Balance  of  momentum  in  the  e2  and  e2  directions  implies  that 

<«■> 


r  _  ^32  , 

C32  —  -j-  +  P  32- 


<221>  ST 

or  equivalently  that 


[3l  +  Ski 

3d 


(2.22) 


(  .  ($31  +  ^32)/  ,> 

r  =  r0(x3,t) - — - (zi,z2,f), 
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0 


whereas  balance  of  momentum  in  the  e3  direction  yields 


(2.23) 


du 

P°W 


dS3 1 
dx\ 


dS32 

dx3 


dlTg 

dx3' 


Here,  p0  is  the  constant  mass  density  of  the  material.  Since  dr0/dx3  depends  on  x3 
and  I],  whereas  all  quantities  on  the  left-hand  side  of  (2.23)  depend  only  on  X[,  x3, 
and  1 1,  we  conclude  that  for  antiplane  shearing  flows  diro/dx3  is  independent  of  x3 
In  what  follows  we  shall  assume  this  quantity  is  zero. 

We  now  turn  our  attention  to  “yield  condition"  and  the  flow  rule  for  the  plastic 
strain  tensor  V  of  (2. 8)2  We  assume  that  yield  is  determined  by  whether  the  scalar 
,S'j,  ■+•  .832  exceeds  a  threshold  ,8^  or  not.  This  assumption  relies  on  the  special  form 
of  5  (see  (2.13))  and  is  equivalent  to  a  yield  criteria  determined  by  the  norm  of 
where 


(2.24) 


2  :  =  SijSij  =  2 (Si  +  S|2)  +  3^(55,  + 


or  one  based  on  the  maximum  shear  stress 


(2.25) 


S: 


def 

:  =  max 


||5e-(5ee)e||2. 


In  the  sequel  we  let  H  denote  the  Heaviside  function 
(2.26)  H(x)  :d=f  {  J’ 


x  <  0 
x  >  0 


and  define  xi>\  and  by 

(2.27) 

and 

(2-28) 


-\L 


H(x  -  5y)  dx 


H(x  —  S,, )  dx. 


where  Sv  >  0  is  the  “yield  stress.” 

We  shall  confine  our  attention  to  the  Coulomb  type  sliding  law 


(2.29) 

dp3\ 

I  dv>\ 

■$31 

dt\ 

fJTg  ds3 1 

“  1JT0 

and 

(2.30) 

dp33 

1  dv<\ 

.833 

dt  1 

>dTg  d.832 

~  Wo 

H ('83i  4-  .832  —  ), 


though  much  of  what  we  say  applies  equally  well  to  the  flow  rule 

SyS3\ 


(2.31) 


dp3\  _  s,j  di>‘ 
dt\ 


JTo  dS3l  jJT0  v/.b'5i  +  .852 


=j=H  f  a/^Si  +  >5'3j  iy 
■831  8 
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and 

(•2..T2) 


dp32  _  Sy  dv? 
&t  1  JTq  di>32 


JTo  \/  “  31  +  -833 


The  constant  3  is  the  shear  modulus  in  (2.12).  Sv  is  the  yield  stiess,  and  T0  >  0  is 
a  fixed  relaxation  time.  The  flow  rule  is  defined  for  Sjj  4-  .Sj2  #  Sy  and  the  problem 
remains  to  define  it  on  the  yield  surface. 

We  first  note  that  if  +Si2  *  5J.  we  can  combine  (2. 17)-(‘2.'20)  and  (‘2.29)  and 
(2.30)  to  obtain  the  following  system  for  631 .  532.  and  u: 


(2.33) 


(2.34) 
and 

(2.35) 


1  OS31 

du 

■"•531  H 1 .731  4-  N32  — 

sV 

3  dt\ 

dx  1 

JTq 

1  $632 

du  _ 

—S3  2  //  ( •<53t  4-  i'5  0  — 

■n) 

3  dtx 

cJxe 

JTq 

du  dS3 1  dS37  _ 

dt  1  dx\  dx7 


Equations  (2.33)  and  (2.34)  imply  that  for  4-  5^  ^  Sj. 


(2.36) 


+ = 23  (531  ik + S^ft 
2 


—  jr('S'31  +  ^32)^(^31  +  5j2  —  Sy), 


and  (2.36),  together  with  the  results  of  [3].  [4],  [5],  motivates  our  extension  of  the  flow 
rule  on  the  yield  surface  5^  +  533  =  Sy .  We  extend  (2.29)  and  (2.30)  to  the  yield 
surface  4-  S3 2  =  Sy  by 


(2-37) 

where 

(2.38) 


dp3i  «s31 


dt, 


JTn 


and 


dp 32  _  O.S32 

dt  i  JTq 


du 


du 


SI 


,  -  -  UU  .  UU  -jv 

1  if  2>3!  +  532  =  by  and  b3i  —  +  ^32—  >  — . 


Oxr 


JTo 


■JTq ( S3 1  j~  +  b'32  ) / Sy  if  63,  4-  .533  —  S~  and 


Of  =  < 


du 


Sr, 


_  Ou  \j  t*  *v 

0  <  2)3!  a - h  ^32  7T -  <  ~p=~  />- 

OXi  0X3  J7q  > 


du 


du 


n  n  n  t* 

0  if  63,  +  .$2 2  =  5;  and  631  - —  4  532- —  <  0e 

OX\  OX o 


4  The  relations  (2.37)  and  ( 2.3d  |  transform  in  a  frame  indifferent  fashion. 
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In  the  sequel  we  shall  coniine  our  attention  to  the  extended  flow  rule  (2. 29),  (2.1(0). 
(2.37)  and  (2.38).  The  relevant  equations  are 


(2.39) 

1  ds3i 

du 

<*•§31 

d  dti 

dxi 

dT0  ’ 

(2.40) 

1  $632 

du 

<*•’532 

d  dti 

dx2 

dT0  ’ 

(2.41) 

du 

dS3\ 

dS32  „ 

PodU  ~ 

dxx 

-  =  0 

OX  2 

where  now 

(2.42) 

1 

if  Si 

+  sh  >  S'- 

1 

if  Si, 

+si  =  s:- 

7 

and 

§(*>: 
y  x 

a  =  < 


PTq  ( r.  du  ,  c  du 
~pr  I  >831  r  032  "X — 

s ;  v  <?*i  <7*2 


°-^(531St  +  S32^ 


if  5|s  +  530  =  Sy  and 


du  \ 

£ 


0  if  Si  +  S32  =  S-  and  ~  (ssi  ^  +  S32j^-  J  <  0  av 


*i 


du  \ 


i  0  if5|1+5|2  <5y2, 


and  these  are  solved  together  with  appropriate  initial  and  boundary  conditions.  Hav¬ 
ing  solved  the  above  system  for  531,  S32,  and  u  we  recover  the  deformation  gradients 
Fu  and  Fz2  by  solving 


(2.43) 


dFzi  du  dFi2  du 

— - - —  —  0  and  — - - — 

Ot  1  UX\  (jt\  ox 


together  with  appropriate  initial  conditions.  The  plastic  strains  p31  and  p32  are  then 
given  by 


A. 

? 

A 

? 


A 

7 


(2.44)  P3,  =  Fu  -  ^7-  and  p3l  =  F32  - 

d  d 

These  equations  should  be  contrasted  with  what  obtains  in  the  more  commonly  stud¬ 
ied  theory  of  rate  independent  elastic-perfectlv  plastic  materials.  In  that  theory  (2.37), 
(2.39)— (2.4 1 ),  (2.43)  and  (2.44)  still  hold  but  a  is  given  by 
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(2.45) 

JTo  (  du  du  \  2  2  2  j 

V  ^7  >!Sn*i„  =  is  .ad 

a  =  <  fL(4i,K7+*!^:)' 

„  ,•>  ,■>  .  5Tn  I  ,  du  ,  du  V 

0  lf.b3t+.b3-2  =  .by  and  _(b31— +  532— )<0. 

.  0  if  +  5|2  <  dy . 

The  unboundedness  of  «  on  the  yield  surface  Sj,  +  532  =  5y  presents  difficulties 
not  encountered  in  our  model.  In  particular,  across  nonstationary  shocks  where  F3I . 
T32 ,  u,  53i,  and  532  experience  jump  discontinuities,  we  must  admit  jumps  in  the 
plastic  strains  p3J  and  p3 2.  The  reason  for  this  is  that  in  the  classical  rate  indepen¬ 
dent  theory — o  as  in  (2.45) — we  must  allow  “dirac”  type  singularities  in  the  terms 
a53i //JTo  and  aSzi/JTo  and  therefore,  we  cannot  conclude  that 

;2.46)  cn![p3i]  =  cn2[p32]  =  0. 

Here,  c  is  the  normal  velocity  of  the  shock  wave  and  n=(ri!,  n3)  is  the  unit  normal  to 
the  shock.  In  our  model  a  is  bounded,  no  “dirac”  type  singularities  arise  in  the  terms 
aSji/fiTo  and  ctSzi/QTq,  and  thus  (2.46)  holds.  This  implies  that  with  our  model  all 
nonstationary  shocks  satisfy  c2  =  1;  that  is,  they  propagate  with  the  speed  of  elastic 
signals.  With  our  model,  the  only  surfaces  across  which  the  plastic  strains  can  jump 
are  stationary,  i.e.  c  =  0.  Such  jumos  are  also  allowed  in  the  classical  theory. 

We  conclude  this  section  by  wnung  down  a  dimensionless  version  (2.39)— (2.44) . 
We  let 


Fill  it-.  -  t — —  J 


PO  *1 

0  To  ’ y  ~ 


P0  £2  .  _  £l_ 

J  T0'  ~  T0 


5U,r31 


and  =  ~ 
7  0 


and  observe  that  (2.39)— (2.42)  transform  to 


(2.48) 


dr3i  dv 
i ft  dx 


(2.49) 


dr32  dv 
dt  dy 


-or32, 


(2.50) 


dv  0r3 


where 


(2.51) 
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1  if  r5l  +  r32  >  Ty  S? 

1  .fril+ri2  =  r;  a„d  + 


f 

dv 

dv 

\ 

(  t31 

^r  +  r32^7 

) 

(  dv 

dv 

1  r2  1 
y 

vr31^i 

+  T32 

dx2 

if  rli 

+  r32  = 

*> 

ry 

and 

if  r3i 

+ 

A 

*> 

The  transformed  versions  of  (2.43)  and  (2.44)  are 

5f3i  du  <9F32  5v_n 

(2-52)  -sr-a;-0  “<i  —  -^  =  0 

and 

(2.53)  P3i  =  F31  -  r3i  and  p32  =  F32  -  r32. 

3.  Uniqueness  results.  Our  task  in  this  section  is  to  establish  the  following 
Theorem  3.1.  Let  fi  be  an  open  domain  in  R2  with  smooth  boundary  dQ..  Then, 
there  is  at  most  one  piecewise  smooth,5  Lfac(Q)  solution  (r3t,r32,  v)  to  (2.48)-(2.5I) 
satisfying 

(3-1)  tlim(r31(r32lw)(r,y,0  =  (r|lt  r&,  t>°)(x,  y), 

(3.2)  (  ljm  (nr-rs!  +  n2r32)(x,  y,  f)  =  fi(z,y,t), 

(r,y)6  H;(x,v)— an, 

(3.3)  Um  a.  *'(*.y.f)  =  /2(aMf.O- 

Here  dQ  =  <9Qt  Uc>n2,  <9Qi  ndfi2  is  at  worst  a  finite  collection  of  points,  n-=(ni,»2) 
is  the  unit  exterior  normal  to  6Qi,  and  the  (functions  figure  smooth  functions  uj 
x  [0, 00)). 

Proof.  We  first  note  that  if  (rj,,  r32,  vb)  and  (r-J,,  r£2,va)  are  two  solutions  to 
(2.48)— (2.5 1) ,  then  their  differences  satisfy 

(3-4)  ~{r 3\  -  Tj,)  —  _  „«)  =  -(^Jt  -  ft-rj,), 

(3.5)  ^(r32  ~  r32)  —  —  V3)  =  — (d4r32  —  d3  r32), 


|(V‘  -  ^  ^  -  T*>  -  ^  -  r3°2)  =  0. 


5  This  formulation  admits  shocks  which  propagate  with  nonual  velocity  c  satisfying  c2  as  1. 
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Here,  d4  and  da  represent  the  bounded  function  d  defined  in  (2.51)  evaluated  at 
(r3i>r32’v4)  aQd  (r3i<  r32' t,a)>  respectively.  The  last  three  identities  imply  that 

2dt  ~  r^2  +  ^ T 22  ~  r^2  +  -  va)i 

(3")  ~  Kr*1  ~  r^vi  ~  ~  ^  ((r32  -  r3a2)(^  ~  »“)] 

=  ~  [(r31  “  T3i)((^[4t31  —  ®ar3l)  +  (r32  ~  r32  ~  d1*^)]  . 

We  now  claim  that 

(3-8)  p  :d=f  (r3\  -  r3a1)(d‘r341  -  d‘r3a  )  +  (r»,  -  r3°2)(d4r42  -  dV3J2) 

is  nonnegative.  In  verifying  this  assertion  there  is  no  loss  in  generality  in  assuming 
that 

(3.9)  0  <  da  <  d4  <  1. 

We  first  note  that  p  may  be  rewritten  as 


(3.10) 


P  —  d°  [(r3i  —  T~3l)2  +  (r32  -  r32)2] 

+(d4  -  d°)  [(r34j)2  -  r31r3l  +  (r42)2  -  r3a2r42]  . 


If  d“  =  0,  then  (r31)2  +  (r3a2)2  <  r2  and  r3a  r3\  +  r3a2r342  <  ry  ^/(r^)2  +  (r362)2  and, 
therefore,  (3.10)  implies  that 

(3.11)  P  >  d^r4,)2*  (r42)2  (V^)2  +  (r42)2  -  ry)  . 

If  d4  =  0,  then  (3.10)  implies  that  p  —  0,  whereas  if  0  <  d4  <  1,  (2.51)  implies 
that  \/(t 3i  )*  +  (r32)‘  ry>  and  (3.11)  then  yields  p  >  0.  We  now  turn  to  the  case 
where  0  <  d°  <  d4  <  1.  If  d4  =  da,  the  nonnegativity  of  p  follows  from  (3.10), 
and  thus  to  complete  the  verification  that  p  >  0  it  suffices  to  consider  the  case  where 
0  <  da  <  d4  <  1.  Here  we  know  that  (T3l)2+(r£2)2  =  r2  and  ('3i)2+(r|2)2  >  r2.  The 
former  identity,  along  with  (3.10)  and  +  r£2r£,  <  ry  ^/(r^,)2  +  (r3J2)2,  implies 

that 


(3.12) 


p  >  da  [(r31  —  r31)‘  +  (r32  -  r|2)'] 

+(d‘  -  dV^4,)2  +  (r342)2(t/(r34l)2  +  (r42)2  -  ry), 


and  (3.12),  0  <  da  <  d4  <  1,  and  (r3l)2  +  (rj2)2  >  r2  complete  the  proof  of  the 
assertion  that  p  is  nonnegative. 

For  any  (i0,  y0)  6  1R2,  r0  >  0,  T  >  0,  and  0  <  1  <  T  we  let 

(3.13)  C'(x0,yo,r0,0  :=f  ((x,  y)|(z  -  x0)2  (y  -  yo)~  <  {ro  +  T  -  l)2}. 

The  identity  (3.7)  implies  that  if  (rj, ,  r^2,  u4)  and  (r31 ,  t£2,  va )  are  two  solutions 
of  (2. 48 )-(2. 51)  taking  on  the  same  data  (3. 1)— (3.3),  then 


(3.14) 


'C(r0,yo,r0,T)nn 


((r3l  —  r31  )2  +  (r32  —  r32)2  +  —  u<>)2)  dx  dy 


^  x 


u  \oV' 

4> 


ANTIPLANE  SHEARING  MOTIONS  OF  A  VISCO-PLASTIC  SOLID  1 

) 

+  [  (  f  .  pix,y,t)dx dy]  dt 

Jo  \yC(xo,»o,fo.T)nn  / 

+  [  (  [  ^  /i((r31  -  r3l)2  +  -  r3a2)J  +  (®*  -  t'a)2)  _ 

Jo  \JaC(x<,,ya.rt.T)nn  \*  ) 

1  (ro  +  T-t)  +  (ro  +  T-t)  j)  ) 


(3.15)  ■  ^C(x0,l/o.r0,0  =  {(x,y)j(x  -  x0f  +  (p  -  yo)J  =  (r0  +  T-I)2}. 

The  vector  ((x  —  x0)/(r0  +  T  —  t),  (y  —  yo)/(ro  +  T  —  t))  is  the  unit  exterior  normal 
to  dC(x0,  y0,  r0,  t),  and  da  is  arc  length  along  dC(iQ,  yo,  ro,  t).  Since 


_(v4  -  v°)  ( zisMlziLI  + 

1  'V  (ro+T-t) 


(y  -  yo)(r32  —  r32)\ 

(ro  +  T-t)  J 


(3.16) 


>  -|v‘  -  va\sj(r^  -  r&)2  +  (r3*2  -  r32 ) 2 

>  -^((r3\  ~  T3ai)2  +  (r|2  -  r32)2  +  (v*  -  i  i-‘ 


and  since  p  >  0,  we  see  that  all  three  integrals  in  (3.15)  are  nonnegative  and  their 
sum  is  zero.  From  this  we  obtain 


(3.!7)  J 


((r3*i  ~  *3i)2  +  (*•»  -  r|2)2  +  (t/4  -  va)2)dxdy  =  0. 


*/  C(ro,yo»^o.T)nfl 

which  is  the  desired  uniqueness  result. 

4.  A  signalling  problem.  In  this  section  we  consider  an  elementary  one-dimen¬ 
sional  signalling  problem  for  the  normalized  system  (2.48)-(2.53).  The  solution  is  of 
the  form 

(4.1)  (r-3'.,r32,u)  =  (r(x.t),0,  v(x,t)),  0  <  x  <  oo. 


where  r  and  v  satisfy 


dr  dv 
dt~fa~  ~aT' 


0  <  x  <  oo, 


^-^1  =  0 
dt  dx  ' 


0  <  x  <  oo, 


'l  if  r2  >  r2  ..-v 

1  if  r2  =  r-  and  >  1, 

A  =  £g£,  if  r2  =  r2  and  0*<  $£  <  1, 

O’  if  r2  =  r2  and  <  0, 


if  r2  <  r2, 
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and  the  initial  and  boundary  conditions 

(4.5)  (r,  v)(x,  0)  =  (0, 0),  0  <  x  <  co 

and 

(4.6)  u(0 ,t)  =  -r0,  where  r0  >  ry. 

We  note  that  the  results  of  the  previous  section  guarantee  there  is  at  most  one  solution 
to  the  above  problem. 

In  the  region  0  <  t  <  x,  we  have  ( r ,  v)  =  (0,0).  Moreover,  r  +  v  is  continuous 
across  the  curve  t  =  x  and  thus  satisfies  r_(f,f)  +  v-(t,t)  =  0.  The  difficult  part  of 
the  problem  is  to  show  there  is  a  curve  t  —  J{x),  Q  <  x  <  x#,  with  -1  <  dj/dx  <  0, 
such  that  in  the  region  x  <  t  <  J(x)  with  0  <  x  <  x#,  r  and  v  satisfy 


(4.7) 

T>  Ty, 

(4.8) 

dr 

ai ' 

dv 

dz 

,  dv 

-t  and  — 

at 

the  boundary  condition  (4.6)  and  r_(M)  +  «-(M)  =  0.  On  the  curve  t  =  J{x) 

we  have  lim<_o+  r(x,  J(x)  -  e)  =  ry  and  v(x)  :=r  Iimf_0+  u(x,  ^(x)  -  «)  satisfies 
0  <  dv/dx  <  Ty.  In  the  region  J(x)  <  t  and  0  <  x  <  x#  we  have  r(x,f)  =  rtJ  and 
v(x,t)  =  D(x),  whereas  in  x#  <  x  <  t,  r  =  ry  and  u(x,  f)  =  v(x#)  =  -ry  (see  Fig,  1). 

The  existence  of  a  curve  t  =  J(x)  with  the  desired  properties  may  be  established 
by  converting  the  system  (4.6),  (4.8),  and  r_(t,t)  +  t>_(<,f)  =  0  to  integral  equations 
for  r  and  v  in  x  <  t,  verifying  that  for  0  <  t  -  x  <  1  the  stress  satisfies  r  >  ry, 
and  finally  by  obtaining  qualitative  information  on  the  level  line  t  —  J(x)  defined  by 
lim,_0+  *•(*,  J{x)  -  «)  =  rv . 
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Rather  than  perusing  that  approach  we  shall  show  how  to  obtain  simple  approx¬ 
imate  solutions  satisfying  (4.6)-(4. 3)  and  r_(t,i)  -f  v_(t,  t)  ss  0  as  well  as  approxima¬ 
tions  to  the  level  line  t  =  J{ x). 

We  note  that  for  each  integer  IV  >  1  the  system  (4.8)  has  solutions 


(4.9) 

and 

(4.10) 6 


VjV  =  -TO  +  A(k)(t)x 


*=1 


■MO  + 

T»=<— - 2k 


where  the  coefficients  satisfy 

(4.1 1)  Aq  +  Aq  =  Ai, 


(4.12)  Afc+A*  =‘2k('2k  +  l)\k+l,  1  <  it  <  iV  -  1, 

and 

(4.13)  V+Av=0. 

These  solutions  satisfy  the  boundary  condition  u.v(0+,t)  =  — r0  and  have  2.V  +  1  free 
parameters  which  are  determined  by  insisting  that  the  equation 

(4.14)  ry{t,t)  +  vy{t,t)  =  0 

is  satisfied  to  0{t‘2N)  as  t  — •  0+.  The  approximate  curve  11  =  Jy(x)  is  subsequently 
determined  by  solving  Ty(x,  JN(x))  =  rv.  An  easy  calculation  shows  that  Jy(x)  - 
O({r0  -  Ty)/T0)  and  djy/dx  <  0  which  guarantees  that  the  number  defined  by 
Jn(x#)  =  x#  is  O((r0  -  Ty)/r0),  and  thus  on  the  boundary  x  =  t,  ry(t,  t)  +  vy(t.  t)  is 
at  worst  0((r<j  —  ry)/ro)2/v+l  for  0  <  t  <  x W'e  continue  the  approximate  solutions 
to  the  rest  of  the  region  described  by  Fig.  1  via  the  extensions  procedure  used  for  the 
exact  solution;  that  is,  for  0  <  t  <  x, 

(4.15)  (riV,u,v)  =  (0,0)  /V 

7 

for  Js(x)  <  t  and  0  <  x  <  x# 

(4.16)  vN(x,t)  =  Vl\(x,Jy(x))  and  r,\{x.t)  =  ry, 
and  for  x#  <  x  <  t, 

(4.17)  Vff(x.t)  =  vlV(x#,  J,\(x#))  and  rN(x,t)  =  ry. 

We  are  then  guaranteed  that  the  error  made  in  failing  to  meet  the  boundary  condition 
ry(t,  t)  +  u,v(M)  =  0  is  at  worst  O((r0  —  ry)/r0)2y+l  for  all  t  >  0.  We  shall  present 
the  details  of  this  procedure  for  the  case  N  =  1. 

In  this  case, 

(4.18)  ui  =  —  To  +  (Ai  o  4-  Aj  ie-1)! 

Here  •  denotes  differentiation  with  respect  to  t. 


a 
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and 

(4.19)  r,  =  (Al>0  +  Altlte-  +  A0lle-)  -  , 

and  the  insistence  that  tq(t,t)  =  0(t3)  as  t  — *  0+  implies  that 

Ai.o  ~r  -'a.i  —  "o 

(4.20)  Alro  -  A0,i  +  2Ai  i  —  0 

Aq,i  —  5Ai  i  =  0 


and  hence  that 


(4.21) 
and 

(4.22) 


t'l  =  T0 


■1  + 


(3  4-  5e-t)x 


.  a  (,+l.-.+ 


The  approximate  curve  t  =  Ji(x)  is  obtained  by  solving  ri(x,  J\(x))  =  ry  or  equiva¬ 
lently  the  equation 

(4.23)  ^3  +  Jxe~*'  +  5e-^‘  - 

The  fact  that  0  <  ry/r0  <  l  guarantees  the  unique  solvability  of  this  equation  for 
0<r<l  and  that  ,7i(0)  =  0(2((ro  —  ry)/ro)).  A  quick  calculation  also  shows  that 


(4.24) 


dJi 

dx 


— 2x 

(8  +  '2Ji  ~  x-) 


<  0. 


The  number  xj^,  where  Jfi(x^)  =  x ^  satisfies 


(4.25) 


3  +  x 


1 


+  5e_r* 


8tv 

To 


and  for  0  <  r0  -  ry  small  enough  we  are  guaranteed  that  x^  =  O((r0  -  Ty)/r0).  This 
estimate,  when  combined  with  (4.24),  implies  that  -1  <  dJxfdx  for  0  <  x  <  x^. 

Our  final  task  is  to  show  that  the  function 


(4.26) 


-  /  \ 

V\(X)  ■-  To 


-1  + 


(3  +  5e-*(e>)*\ 

5  ) 


satisfies 


(4.27) 


0  < 


dv  1 
dx 


(z)  <  Ty , 


0  <  x  <  x^. 


The  defining  relation  (4.26)  implies  that 


(4.28) 


dv  1  /3  +  5e""r,<x) 

— W  =  - 5 - 


5r0e-'7lfr> 

8 


*j;(*). 
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and  this  relationship,  when  combined  with  (4.23)  and  (4.24),  implies  that 


/..on.  _  ,  Toe-J'^x2 -2Ji)($-r2Jx  -  r2)  +  20r2) 

(4..9)  *■(*>-•.+  - - 

The  fact  that  J\{x)  >  for  0  <  x  <  =  O((r0  —  ry)/r0)  implies  that  the  second 

term  in  (4.29)  is  negative  and  this  provides  the  desired  upper  bound  for  di\/dx.  The 
desired  lower  bound  is  an  immediate  consequence  of  (4.2 S)  and  tne  bounds  for  djxjdx. 

We  conclude  this  section  by  contrasting  the  above  solution  with  what  obtains  if 
we  replace  our  flow  rule — o  given  by  (4.4) — with  the  one  generated  by  (2.31)  and 
(2.32)  and  also  by  the  flow  rule  associated  with  a  rate  independent  elastic-perfectly 
piastic  material.  In  the  former  case,  (4.2)  is  replaced  by 


(4.30) 


dr  dv 

dt~  dx~  ~aTy' 


and  (4.4)  is  unchanged. 

In  the  region  0  <  t  <  x  we  have  (r,  v)  =  (0,0),  and  r  +  v  is  continuous  across 
x  —  t.  For  0  <  x  <  t  <  2(tq  —  ry)/ry  we  have 


(4.31)  v  =  -r0  +  ^  and  r  =  r0  - 
for  0  <  x  <  2(r0  -  ry)/ry  and  t  >  2(r0  -  ry)/ry  we  have 

(4.32)  v  —  —tq  +  and  r  =  ry, 
and  finally  for  2(r0  -  ry)/ry  <  i  <  f  we  have 


(4.33) 


v  =  —  ry  and  r  =  ry. 


With  this  flow  rule  the  curve  t  =  J( •)  is  the  constant  function  J(x)  =  2 (r0  - 
ry)/ry,  0  <  x  <  2(r0  — ry)/ry.  Equations  (2.52)  and  (2.53),  the  initial  conditions 
(fr3i,P3i)(x,0)  =  (0,0)  for  x  >  0,  and  (4.31 )— (4.33)  allow  us  to  determine  {Fzi,P3\)- 
The  result  is 


(4-34)(F3i,P3.)=  < 


(0,0),  0<f<x, 

(ro  +  ry(|  ~  x),  ry(<  —  x)),  0  <  X  <  t  < 

(to  +  Ty(|  -  x),r01  -  Ty  +  ry(|  -  r)),  —lltl  <  t  and 

o  <  x  < 

(r,.0  ),Hil^J<x<I. 


It  is  worth  noting  that  the  above  solution  is  unique.  This  can  be  established  using 
the  arguments  of  §3  directly  on  the  system  (4.30)  and  (4.3)-(4.6). 

We  now  examine  the  signaling  problem  for  a  rate  independent  elastic-perfectly 
plastic  material.  Equations  (4.1 )— (4.3)  and  (4.5)  and  (4.6)  still  hold,  except  now  a  is 
given  by 


r  9v 

{  ^ 


if  r2  <  r2 

if  r2  =  r2  and  <  0, 

ifr2  =  r2  and  0<^||. 


(4.35) 
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dFzx  dv  dpil 

dt  dx  ’  dt 


=  ar,  and  F31  =  r  + 


and  these  satisfy  the  initial  conditions 

(4.37)  (F3,,p3i)(x,0)  =  (0,0),r>0. 

We  seek  solutions  with  structure  similar  to  that  obtained  for  the  previous  two  models. 
Specifically,  a  shock  curve  t  =  t(x)  such  that  in  the  region  0  <  i  <  t(x), 

(4.38)  (F31,P3i,r)U)  =  (0,0,0,0), 
and  in  the  region  t  >  i(x)  the  shear  stress  r  is  at  yield,  i.e., 

(4.39)  r(x,  f)  =  ry,  t(x)<t. 

We  interpret  (4.3)  and  (4.36X.3S  conservation  laws,  and  this,  together  with  (4.38)  and 

(4.39) ,  implies  that  on  l  =  t(x), 


(4.40) 
and 

(4.41) 


t/“(x,f(x))  +  ry^~  =  0 


Here,  ( v~ ,  F3^)(x,  l(x))  =  Iime_0+(i>,  F3i)(x  —  (,  t(x)).  The  identity  (4.39)  also  implies 
that  in  t  >  t(x)  the  velocity  v  is  a  function  of  x  only.  Near  i  =  0«e  choose 

(4.42)  v(x,t)  =  -r0  +  Ax,  A  >  0. 

With  this  choice  we  obtain 

(4.43)  P31  =  A(f  -f(x))  +  p_(x) 
and 

(4.44)  F3l  =  tv  +  \(t  -  t(x))  + p-{x). 

Equation  (4.40),  together  with  1(0)  =  0,  then  yields 

Tq  -  (r0  -  Ax)2 


(4.45) 


t(x) 


2Ar„ 


and  (4. 41), (4. 44),  and  (4.45)  imply  that 


/*; 


(4.46) 


p_(x)  = 


(r0  -  Ax2)  -  r2 


We  now  let 
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and  note  that 

(4.48)  p_(z)  >  0,  0  <  z  <  *#. 


(4.49) 


P-(*#)  = 


and 


(4.50)  g|(*#)  =  I- 

In  the  region  (rj  —  (ro  —  Ar)2)/2Ary  <  t  and  0  <  x  <  z#  =  (ro  —  ry)/A  our  solution 
is  given  by 


(4.51) 

(4.52) 


(4.53) 


v  =  —  ro  +-  Az, 


(4.54)  r  = 
The  shock  curve  is  continued  to  x  >  x#  by 

r 2  -  r 2 

(4.55) 


and  in  the  region  ((r|-r2)/2Ary)-)-(i-((ro  -ry)/A))  <  t  and  (r0-ry)/A  =  r#  <  x, 


(4.56) 


^31  =  ry ,  P3i=0,  v  =  -ry,  and  r  =  ry 


The  line  x  —  x#  —  (r0  —  ry)/A  is  a  stationary  contact  discontinuity  and  across  it  P31 
jumps  while  the  other  fields  are  continuous.  The  interesting  fact  about  the  signaling 
problem  for  this  model  is  the  lack  of  unicity  of  solutions;  we  have  a  compatible  solution 
for  every  A  >  0.  This  observation  points  out  one  of  the  weaknesses  of  the  classical 
model. 

5.  Computational  experiments.  In  this  section  we  present  the  results  of  a 
computational  experiment  performed  on  the  normalized  system  (2.48)-(2.52)  when 
the  pressure  gradient  is  zero.  The  results  reported  deal  with  a  two-dimensional  gen¬ 
eralization  of  the  signalling  problem  of  the  previous  section. 

The  experiment  deals  with  the  system  (2.48)— (2.51)  solved  in  the  region  r  >  0 
and  t/2  <  9  <  2t,  where  r  =  \/i2  +  y- .  At  time  t  =  0  we  assume  that 


(51)  foi,  nj2. »)  =  (0,0,0) 

for  r  >  0  and  x/2  <  9  <  2t,  and  for  t  >  0  we  assume  that 

v  fr,  ^  j  =  t>(r,2x~)  =  r0> 


(5.2) 


r  >  0, 
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Fig.  2. 


where  rQ  >  ry,  and  again  ry  >  0  is  the  yield  stress. 

The  elastic  version  of  this  problem,  namely,  the  system 


(5.3) 

di-31 

dt 

(5.4) 

<9r32 

dt 

(5.5) 

dv 

OT31 

Hi 

dx 

dx 

dv 

dy 


=  0, 


T32 


dy 


=  0, 


together  with  (5.1)  and  (5.2),  was  considered  by  Keller  and  Blank  [13].  They  obtained 
exact  solutions  to  this  and  a  number  of  other  problems  with  self  similar  structure. 
Relevant  to  us  here  is  the  singular  nature  of  r|j  +  r|2  as  r  — *  0+.  Their  results 
demonstrate  that 

(5-6)  Tjj  -f  r32  =  O  ,  r  — *•  0+. 

This  singular  behavior  also  obtains  for  the  plastic  flow  problem  and  forces  us  to 
treat  the  boundary  conditions  in  our  numerical  simulation  carefully.  Our  integration 
scheme  for  (2.48)-(2.51)  is  based  on  a  symmetrized  operator  splitting  algorithm  for 
the  governing  differential  equations.  At  time  t  =  nh,  n  —  0,  l,  2, . . our  approximate 
solution  consists  of  lattice  data 


(5.7)  (r31l  r32,  v)(*fc,m)  =  (1*3, ,  r32,  v)  . 12/t|  nh^  . 

For  the  problem  under  consideration  the  boundaries  are  not  part  of  the  computational 
lattice  but  are  offset  from  it  by  a  distance  of  A/2.  The  computational  lattice  is 


(5.8)  S  =  {(k,  m)  j  k  <  0  and  m  =  0,±1,  ±2, . . .}  U  {(*,  m)  |  k  >  1  and  m  <  0}  . 


To  update  the  data  (5.7)  we  successively  solve 


di-31 

m 


dv 

dx 


0, 


^r32 

dt 


—  0,  and 


dv 

Hi 


dr3\ 

dx 


=  0,  0  <  t  <  A, 


(5.9) 
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d 


.  i,r.- 

r  ,w 


(5.10) 

<^31  0f32 

dt  ’  dt 

dv 

— —  =  0,  and 
dy 

dv 

dr32 

dy 

and 

(5.11) 

dr3i 

dt  =  aT3!l 

dr32 

nr  = 

and 

II 

=  0,  0  <  t  <  h, 


=  0,  0  <  t  <  h, 


where  of  course  a  is  defined  in  (2.51).  For  (5.9)  we  use  the  approximate  solution 
defined  by  (5.7)  as  initial  data  and  let  r32>  t'Mlt.m)  denote  the  value  of  this 
solution  at  t  ss  h  on  the  lattice  $.  We  then  solve  (5.10)  using  the  (rj,  r^2,  td  ),  *  ml 
as  initial  data  and  let  (r|()  r|2.  irW ,m>  denote  value  of  the  solution  at  t  =  h  on  S. 
Finally,  we  solve  (5.11)  with  (rli.fjj,  tr)(jfc  m)  as  initial  data  and  let  ( ■  r32  ■  y3)r*  mj 
denote  the  value  of  this  solution  at  t  —  h  on  5. 

We  then  repeat  the  process  solving  (5.10)  first  with  the  data  (5.7),  and,  we  let 
(r3i  t  r32>  u4)(*,m)  denote  the  lattice  update  at  t  =  h.  We  then  solve  (5.9)  using 
(r3ii  r32’  u4)(t,m)  33  initial  data  and  let  ("It , 7*32,  ws)(i  m)  denote  the  lattice  update. 
Finally  we  solve  (5.11)  with  data  (t^  ,  r|2,  u5)(t  m)  and  let  (fit ,  r|2,  t/6)(t,m)  denote 

the  lattice  update  at  t  =  h.  The  desired  approximate  solution  (731, 732,  t/)^^  is  then 

and  (f|t,f362>  v6)(i,m);  that  is, 


obtained  by  averaging  (f|t ,  f|2,  v3)(t,m)  *uu  v '31X32 


(5.12)  (fai.fsz.w^t!)  =  j(f; 


3  ,  T« 
31  +  r31> 


32 


■3i> 


f3  +  V6)(t,m). 


Of  course,  all  of  the  intermediate  updates  are  solved  subject  to  the  boundary  con¬ 
ditions  of  the  original  problem.  Here  these  boundary  conditions  manifest  them¬ 
selves  as  reflection  conditions  at  those  lattice  points  that  are  a  distance  h/2  away 
form  the  actual  boundary.  Formal  accuracy  could  be  maintained  if  we  used  either 
(r|v,  r32>  v3)(k,m)  or  (r!i>  r32>  for  the  updated  approximate  solution  but  either 

of  these  updates  alone  would,  over  time,  tend  to  introduce  asymmetries  into  the  ap¬ 
proximates  not  present  in  the  actual  solution.  These  asymmetries  are  removed  with 
the  algorithm  employed. 

The  results  of  our  experiment  are  shown  in  Figs.  3-7.  Each  snapshot  shows  two 
different  representations  of  the  velocity  field  and  the  total  shear  stress,  namely  the 
quantity  s/r3i  +  This  simulation  was  run  with  h  =  i/50,  ry  =  1.  and  r0  =  1.3./ 
The  contours  on  the  velocity  plots  are  spaced  0.1  apart  and  run  from  v  =  0  to  v  =  1,3  J- 
The  stressjcontours  run  from  1  to  3.2  in  increments  of  0.2.  In  these  snapshots^gj^ 
jiolTsee  any  of  the  plane  wave  solutions  of  the  previous  section  but  also  the  effect  of  V 
the  corner  singularity  which  are  confined  to  the  region  0  <  r  <  t  and  t/2  <  9  <  2x.  \ 

For  comparison  we  have  run  the  elastic  version  of  this  problem  with  the  same 
boundary  conditions  and  same  values  of  h,  ry,  and  r0.  These  results  are  shown  in 
Figs.  8-12. 

It  should  be  noted  that  for  both  problems  the  velocity  fields  satisfy  the  additional 
condition 


r> 


Iim  v(r,9,t)  =  0, 

r — t*1 


-  <  9  <  2tc 


(5.13) 

and  that  our  numerical  solutions  meet  this  consistency  condition  automatically. 


r- 


IMWIM 


MM 
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velocity 


dmc0.8 

maximum  total  shear  stress4.143 


total  shear  stress 


total  shear  stress 


Fig.  11 


timel 

maximum  total  shear  stress4.404 


total  shear  stress 


total  shear  stress 
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